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Abstract 

In this paper we show that the standard causality condition for atten- 
uated waves, i.e. the Kramers-Kronig relation that relates the attenuation 
law and the phase speed of the wave, is necessary but not sufficient for 
causality of a wave. By causality of a wave we understand the prop- 
erty that its wave front speed is bounded. Although this condition is 
not new, the consequences for wave attenuation have not been analysed 
sufficiently well. We derive the wave equation (for a homogeneous and 
isotropic medium) obeying attenuation and causality and with a gen- 
eralization of the Paley- Wiener-Schwartz Theorem (cf. Theorem 7.4.3. 
in [5]), we perform a causality analysis of waves obeying the frequency 
power attenuation law. Afterwards the causality behaviour of Szabo's 
wave equation (cf. [llj ) and the thermo- viscous wave equation are inves- 
tigated. Finally, we present a generalization of the thermo-viscous wave 
equation that obeys causality and the frequency power law (for powers 
in (1,2] and) for small frequencies, which we propose for Thermoacoustic 
Tomography. 
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1 Introduction 

In physics an attenuated wave is modeled by replacing the real frequency de- 
pendent wave number by a complex frequency dependent wave number. For an 
attenuated spherical wave with origin (x, t) = (0,0), this yields [3"1 fTU], 17] 

1 r e -a(a))|x| g-tw^-J^y) 

(1) Pq(x ' t) = v^J d " ' 

R 

where a = a(u>) > is called the attenuation law, which is assumed to be a 
positive real-valued even function, and c = c(cj) is called the phase speed. The 
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complex wave number is given by k(u>) = + i a(uj). With the notion 

(2) a* (u) := a(oj) — i [ — ^— ] = — i k{u) + i — (c > const.) , 

\C{U) Cq) C 

the standard causality requirement for waves of the form {!]) read as follows 
(cf. [T3] and Section 3.3 in [T]) 

(3) Re(a,(w)) = -lm(H(a*(u))) and Im(a»(w)) = Re(W(a*(w))) . 

These relations are the Kramers- Kronig relations for a* and are equivalent to 
the relation (cf. [T3"] ) 

(4) -y-r --=-W(a(w . 
v ; c{uj) c 

Here Ti. denotes the Hilbert transform. 

In this paper we require that every wave has a finite front speed v a . Of 
course, this is not a new requirement (cf. [2]), but the consequences for wave 
attenuation have not been analysed sufficiently well. If a wave p a of the form (JTJ) 
has a finite front speed v a bounded from above by a constant vb > 0, then the 
distribution 

Ixl i . 



9a(x,t) :=47r|x|p a \^,t+^-J x € WC, t € K 

is causal for any x € R 3 , i.e. its support lies in [0, oo) for any fixed x e R 3 . 
Conversely, if g a (x,-) is causal for any x £ R 3 , then the front speed of the 
wave p a is bounded from above by vb > 0. Let g a {x.,u>) denote the Fourier 
transforn^ of g a (x, t) . Then g a is causal if and only if the Kramers-Kronig 
relations for g a holds (cf. Section 3.3 in pQ), i.e. 

(5) Re(g a ) = -Im(H(g a )) and Im(g a ) = Re(H(g a )) , 

If g a (x, •) is a Schwartz function then these conditions imply the Kramers-Kronig 
relations © for a* . This can be seen as follows: Let x S R 3 be arbitrary but 
fixed and let *t denote the time-convolution. For a spherical wave of the form ([T]) 
with Co := vb we have 

(6) ffa (x,t) = -L^(e-«*^l x l)(i) (xeR 3 , t6R) 

V 2 7T \ ' 

with a* defined as in (|2|). Since the causality of g Q (x, •) implies the causality of 

(7) Vff a (x,i) = -— a*(i) *t3 Q (x,t), 

a*(t) must be causal. But the causality of means that the Kramers-Kronig 
relations j3]) for a*(u;) hold. We note that in general, the causality of d* or Vg a 
does not imply the causality of g a (cf. Section d]). 



1 If / e 5(K), then />) := := ^ / R e^" *f(t)dt and / := ^(Z). This 

convention implies fg = -jL= f * g and fg= -]L- f * g. 
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The phase speeds induced by the frequency power laws a(u>) — ao M 7 with 
7 G (R + \N) U {1} and the Kramers-Kronig relation are derived in [TBI [TBI 
rrSl [T^l [T2] . For example if 7 G (0, 00) and 7 G" N, then the phase speed reads as 
follows 
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c(w) c(lo ) 
which implies 



a tan(| 7 ) flwp- 1 - |"o| 7_1 ) (w G K) 



(8) a*(u/J = + za tan I -7) cj 7 w. 

cos(-2 7j V2 / 

In Section 0] it will be shown that the wave (fTJ) with this a* is causal only if 
ujq = and 7 G (0, 1). This shows that the Kramers-Kronig relation |[3| for a* 
is necessary but not sufficient. Moreover, this shows that the frequency power 
law for 7 > 1 is not an admissible attenuation law. 

A standard calculation shows that the wave ((TJ) with a* defined as in (JSj) 
satisfies the following wave equation: 



(9) VV 



cos(f 7) * cq dt 



p a = -S(x)S(t) , 



where ao > and Dj denotes the Riemann-Liouville fractional derivative de- 
fined by (0 El)) 

(io) r'lD/f/DH^H^/H. 

We note that the causality of the kernel of D] implies the Kramers-Kronig 
relation J3J and vice versa. 

A different wave equation was derived by Szabo in [llj). A comparision of 
the dispersion relations of the ld~ thermo-viscous wave equation and the ld- 
electromagnetic wave equation for conducting media lead Szabo to the following 
dispersion relation k(uj) 2 = ^+i 2 — ao M 7 with 7 > 0, which implies the wave 

c co 

equation: 



For example, if 7 > and 7 G" N, then 

(12) L 7 = - 2 °° f(H<). 

v2tt cos(| 7) c 

Since L 1 is the kernel of the operator cos ^°) CQ , Szabo's equation can 

a 2 2 2 

be obtained by neglecting the term cos 2°'k ^ D t 7 in equation (|9|) . In Section [5] 
we show that Szabo's equation admits causality only if 7 G (0, 1). And only 
in this case Szabo's equation is a small frequency approximation of equation Q. 



This paper is organized as follows: The general properties of attenuated waves 
assumed in this paper are introduced and discussed in Section O With these as- 
sumptions we derive the general wave equation for a homogeneous and isotropic 
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medium obeying attenuation and causality in Section^ The first causality anal- 
ysis of attenuated waves is performed in Section [3] for the case of the frequeny 
power law. Afterwards the causality analysis of Szabo' equation (Section[5]) and 
the thermo-viscous wave equation (Section are perfomcd. Finally, a gener- 
alization of the thermo-viscous wave equation, is derived in Section [7j which 
admits causality and obeys the power frequency law for sufficiently small fre- 
quencies and 7 G (1,2]. 

2 General properties of attenuated waves 

In this section we postulate the basic properties of attenuated waves propagating 
in homogeneous isotropic media and infer the structure of attenuated waves. 
The main goal of this section is to clarify our assumptions on which we base our 
casuality analysis of attenuated waves. 

Assumptions 

Let po denote the solution of the standard wave equation V 2 po — \ ^iSr = — f 

c a at 

with po|t<o = and ^§f\t<o — 0. Let A a denote the map that associates to 
each source term / the corresponding attenuated wave p a . a = means no 
attenuation, i.e. pq = Ao(f). We assume a homogeneous isotropic medium, 
which implies that 

A a commutes with all translations in x and t . 

We also assume that for a ^ 0: 

A a is a linear continuous mapping 
that maps C£°(IR 4 ) into S(R, C°°(R 3 )) , 

and that 

Aa{S{x)S{t)) e 5(M,C(R 3 )). 

Here S denotes the space of Schwartz functions, 6{x) denotes the 3D— delta dis- 
tribution and S(t) denotes the ID— delta distribution. The last two assumptions 
take into account that wave attenuation smoothes and decreases the wave. 

Superposition law 

Since the operator A a satisfies the assumptions of Theorem 4.2.1 in [5J, A a is 
a space-time convolution operator with kernel G Q (x, t) :— A a (S(x) S(t)). Since 
5(x) 5(t) is the source term of a spherical wave for a = 0, we interprete G a as 
an attenuated spherical wave and thus 

(13) A*(/)(x,i) =G a (x,t)* x , t /(x I i) 

is nothing else but the superposition law of attenuated waves. In analogy to 
partial differential equations and linear system theory we call G a the Green 
function of the wave model. 



4 



Causality condition 

We require that any attenuated spherical wave G a has a positive finite wave 
front speed v a , which is equivalent to the requirement that 

(14) g a (x,t) :=4 7r|x|G Q (x,i + T(x)) (v a > 0) 

is a causal distribution i.e. supp(g Q (x, ■)) C [0, oo). HereT(x) := Jjj - ^ dr > 
denotes the travel time of G a . Then any attenuated wave A a (f) with com- 
pactly supported source / has finite front speed if and only if the Green function 
G Q (x, •) has finite front speed. We note that property (fl4"|) implies 

supp(G a (x, •)) C [0, oo) for all xel 3 . 

Moreover, we assume that the front speed v a is continuous, which implies T(x) G 
G(R 3 ) and g a 6 £>(R, G(R 3 )). For the standard case were the wave front speed 
is assumed to be constant, this assumption is satisfied. 

General structure of attenuated waves 

The property g a S 5(R, G(R 3 )) implies that 

(15) g a ( X ,t) = ^=F(e-^M) , 

with Re(/3*) := — log (v^7r|ff a |) an d Im(/?*) := — arg(g Q ). Here g Q (x, •) denotes 
the inverse Fourier transform of g a . Since g a is real- valued, it follows that 
Re(/3*) is even and Im(/?*) is odd with respect to oj. (We note that the causality 
of <7a(x, ■) implies that z G R+iM + i— > g Q (x, z) £ C is analytic (cf. Theorem l4.1|) 
and thus — Re(/3*(z)) : R + «R + [— oo, oo) is subharmonic (cf. Example 4.1.10 
in 0).) 

If ^rlm. w)) < - , then /3* can be written as follows 

(16) fl.(|x|,w)=o 5 (w) + j 



where a and c are positive and even with respect to lu and lim^^oo Re(/3*) = oo. 
Moreover, as is a complex valued function such that 

g a (0,t) e5(l). 

By we denote the time-convolution operator with kernel g a (0,t) that maps 
S(M.) into <S(R). If depends only on uj, i.e. /3* = 05, then the attenuated 
wave satisfies the standard wave equation with source term / replaced by S a f. 
We note that g a is causal if and only if Theorem 14. II (cf. Section [4]) is satisfied. 
The standard normal form ([T]) assumes that 

S a = Id, v a = cq = const, a — a(oj) and < a(u>) < 00 , 

which together with the Kramers-Kronig relation |4]) imply that c = c(lu). More- 
over, /3*(|x|,u0 = a*(ui) |x| with a* defined as in @. We note that g a is causal 
if and only if a* satisfies Lemma 14.21 



a(r, lu) 



c(r,w) v a (r) 



dr . 
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Since wave attenuation is an irreversible thermodynamic process, /3* can de- 
pend on time or equivalently on the distance to the origin. We will see that in 
general, wave attenuation depends on its history, even if a depends only on lu. 
In particular this is the case if the support of the kernel of S a is not discrete. 
Moreover, it is not evident that the front speed of an attenuated spherical wave 
is constant. Therefore we see no reason to assume that /3* and the front speed 
v a depend only on u. 



Remark 2.1. If v a (r) < vb for all r > (vg = const.) and g a is a causal 
distribution, then g(x,t) := 47r |x| G a ^x, t + j^j-^ must also be a causal distri- 
bution. In this case a* corresponding to <?(x, t) is given by hi 6}) with v a replaced 

by vb and the Green function reads as follows G Q (x, t) — g a (x,t) * t 4 ^ |xj ^ • 
This fact will be used to prove the non-causality of some wave models. 



3 Wave equation obeying attenuation and causal- 
ity 

Now we derive the wave equation satisfied by the attenuated waves described 
in Section [2] and discuss its Cauchy problem. 

First we derive the wave equation for the Green function G a . The most 
convenient derivation uses the representation of the Green function introduced 
in Remark 1 2. II 



(17) 



G B (x,i) 



G a = g a *t Gb 

*(*-£) 



with 



4 7T Ixl 



and 



V £ 7T \ ' 



Here the constant vb is an upper bound of the front speed of G a and /3* is 
defined as in (flU|) with w Q (x) replaced by vb- We recall that S a is the time 
convolution operator such that S a S(t) = g a (0,t) E S(M). To formulate the 
wave equation we need the time convolution operators : V + (R) — > V + (M.) 
and D*,\tc\ '■ T>' + (W) — > T>' + (M) with causal kernels 

(18) K,(r,t) = -^=^(^(r,uj)\(t) and K^ x] (r,t) = ^M) 

for all r, t > 0, respectively. Here V, (R) denotes the space of causal distribu- 
tions. From (fTT|) . it follows 



K* *t g a 



This together with (fTT|) imply 
1 d 2 G a 



|.t| 



at 2 



and 



D* — ZL 



kl 



G Q -S tt <5(i)<5(x). 
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Due to causality of g a we have 



(19) 
From 



Ga 



t<o 



dt 



dga . - aw. /■■•;%:■: 

9a ' OX 



,V 2 9a e5(M,C(K 3 )) fori = 1,2, 3, 



it follows that 



Dt 



_2_ d_ 

vb dt 



D*-D 



*>\ x \ 



G a 6 5(K,C(R 3 )) 



and that the space-time-convolution operator A a with kernel K\ a maps 5(K, C°° (M 3 )) 
into <S(K, C°°(M 3 )). Since an arbitrary attenuated wave is of the form p a = 
G a *x,t /, the general wave equation reads as follows: 



v v - 4- = ( A « - w) with 



(20) 



v% dt 2 



Pa\t<0 ~ ^ 1 



dp a 

dt 



= 0. 



t<0 



This equation has for every / 6 <S(R, C 
solution 



with compact support a unique 



4 7T |x — X' 



dx' e 5(R,C°°(R 3 )) 



with finite wave front speed. 

If the attenuation law and the phase speed do not depend on the spatial 
position, then the operator does not depend on the spatial position, too, 
and 



(21) 



D*(G a ) * X)t / = £>*(pa) 



holds. Moreover, then £)* ui is the zero operator. In this case we can write the 
wave equation as follows 



(22) 



V 2 p„ 



S a (J), 



with p a | t<0 = and -§f 



t<o 



0. In general the supports of the kernels of D* 



and S'q, are subsets of [0, oo) with positive Lebesgue measure, which means that 
the attenuated wave depends on its history. Since the values of the wave in the 
past are required, the Cauchy problem of wave equation l]22[) is not reasonable. 
In the next theorem we formulate a generalization of the Cauchy problem and 
state its properties for the special case a = a(uj), c — c(co) with as = 0. 
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Proposition 3.1. Let be the time-convolution operator with causal kernel 

K« defined as in h!8\) and let q G <S(R 3 x Rg ), (p :— <z| t=0 and tp := 

Moreover, let G a denote the Green function of wave equation I2ty) with S a = Id 
and wave front speed < vb ■ Then the solution of the generalized Cauchy problem 

2 



vV - 

Pa\t<0 : 



vb at 



Pa=0 



and 



dp a 
dt 



1>, 



t=o 



is given by 
with 



Pa — G a * X) ( / 



V B V B 



Dt + — ^-D*,M m 
vb at w 



q. 



Here Mjjm denotes the multiplication operator, H — H(t) the Heaviside func- 
tion and [,] denotes the commutator, i.e. if A, B are operators then [A,B\ = 
AB-BA. 

Proof. Let p a , G a , f and H (t) be denned as in the Proposition. For convenience 
let A* := D 2 + Then p a := H(t)p a satisfies the following properties: 



d 2 p Q 
dt 2 



= H(t) 



V 2 p a = H(t) V V , 

2 p a 



dt 2 



i>5{t) + v 8\t), 



and 



A*p a = H A* (p a ) + [A* , Ms] Pa , 



since ip = p a \ t=0 and ip = -|» 
1 d 2 p a 



From these properties we infer 



V 2 p« - ^ 



Pa\t<0 



V% dt 2 



= and 



A*p a = —-t- 8(t) — %- S'(t) — [A*, Mff] Pa j 
V B V B 

dp a 

dt 



= 0- 



t<0 



which has the solution p a = G a * x ,t / and thus p a = G a * x ,t / on 
This proves the Proposition. 



□ 



4 Causality analysis of attenuated spherical waves 
obeying the frequency power law 

In this section we show that not every attenuated spherical wave of the form ([1]) 
has a finite front speed although its phase speed is properly related to the 
frequency power law via the Kramers-Kronig relation ([¥]). Other attenuated 
spherical wave models are analysed in Section and [7| 

In the following we use the notions R + — (0, oo), Rj = [0,oo), R~ = (— oo,0) 
and Rq" = (—oo,0]. The next Theorem is a reformulation of Theorem 7.4.3 
in [5] for the case of causal tempered distributions. 
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Theorem 4.1. A distribution u £ <S'(R) is causal, i.e. supp(u) C M.q , if and 
only if 



(Al) 



C, z n u{z) is analytic and 



(A2) 3e > 03C> 03N > OVz £ R + i (-00, -e) : |6(z)| < C(l + \z\) N . 

Applying Theorem 14.11 to attenuated spherical waves of the form (fTJ) yields 
the following Lemma. 

Lemma 4.2. Let p a be defined by (QP with real-valued functions a — a{uj) > 
and c = c(tu), and let p Q (x, ■) £ S'(M) for any x £ R 3 . The wave defined 
by (Qp has a finite wave front speed if and only if a* defined by (0) satisfies the 
following conditions: 

(Bl) a* = z) is analytic onR + iR^ and 



(B2) Be > 03C > 03N > OVz £ R + i(-oo,-e) : 
- fle(a*(-«)) < C + W log(l + |z|) . 

Proof. According to Theorem 14. II a„ fx. •) defined as in (fT4"| is a causal distri- 
bution for any fixed x £ R 3 , if properties (Al) and (A2) hold for u := g Q (x, •) 
for any x £ R 3 . Let x £ R 3 be arbitrary but fixed. According to ^ we have 
5 Q (x, •) = T [er a '^ l x l) = T~ x ( e - Q *(-") l x l) and thus 



(23) 

Let /3i(z) := - 
is analytic on 



(24) 



-Re(a*(-z)) |x| and /3 2 (z) := -Im(a,(-z)) |x|. Since e ~ Q *(~ z )l x l 
R + i Mr , the Cauchy-Riemann equations are satisfied 

'2h. _ 2M cos/ 3 2 = \^i + 2^i 

dx dy dy dx 

_ sin a 2 = _ & + ^a" 

dx dy dy dx 



sin fa , 
cos /?2 



which imply for all z with cos/32(z) ^ and sin ^(z) ^ the equations 

dfc^dfh and 3ft = 

dx dy dy dx 

This means that (3\(z) + i ^(z) is analytic for all z £ R + it" satisfying 
cos #2 (2) 7^ and sin ft (z) 7^ 0. The same equations follow easily if sin ft (z) 7^ 
but cos ft (z) = 0, and for cos ft(z) 7^ but sin ft(z) — and thus a*(— z) is 
analytic on R + «R~. This shows that condition (Bl) is satisfied. Conversely if 
a* (— z) is analytic on R + i R~ then due to the chain rule e~ a '^~ z ^> l x l must be 
analytic on R + iM.~ , since the complex exponential function is analytic on C. 
From condition (A2) in Theorem 14.11 together with ([23l we infer 

e -Re(a.(-»)|*|) <c{l + \z\ 



whereby we can assume without loss of generality C > 1. Let C :— logC > 0. 
Since the real (natural) logarithm function is monotonic increasing, we can apply 
it onto the previous inequality, which yields condition (B2) with C replaced by 
■pr. Conversely, if (B2) holds then condition (A2) holds, too. This concludes 



|x| 

the proof. 



□ 
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For the analysis of the causality properties of the standard models of spher- 
ical attenuated waves we need the following Lemma. 

Lemma 4.3. Let x G K 3 and 7 > be arbitrary but fixed. Moreover, let 
s(7) be the sign o/cos(|7) i/cos(^7) ^ and let s(-f) — 1 i/ cos(| 7) = 0. 
The distribution ioj) 1 ) is causal on K for any 7 > and the distribution 

jF(e~ s ( 7 ) ( -lL,J ) T l x l) is causal if 7 G (0,1] and non-causal if 7 > 1. 

Proof, i) First we prove the causality of J r ((— i w) 7 ). Since w£Ih (— i cj) 7 G C 
is a slowly increasing function, it is a tempered distribution and hence 

i w) 7 )(i) 6 K is also a tempered distribution. Therefore Theorem 14.11 
can be applied. Let 7 € (0, 00) be arbitrary but fixed and let J7 denote the 
complex plane without the non-negative real axis. Then z 7 := e 7 log2 is analytic 
on U (cf.[8]). Therefore (i z) 7 is analytic on R + iIR~ and condition (Al) in 
Theorem I4.1l is satisfied. Moreover, (iz) 7 G C satisfies condition (A2) 

in Theorem 14.11 with C = 1 and TV = 7, which shows that ^ r ((— i aj)~ f )(t) is 
causal on K. 

ii) For the rest of the proof let x G M 3 be arbitrary but fixed. Now we prove 
that V2 7r g a (x, t) := J : [e~ s ^^ luJ ^ ^) is causal for 7 G (0,1] and non-causal 
for 7 G (l,oo). Since (-iw) 7 = |w| 7 [cos( : |7) — isgn(w) sin(|7)], it follows 
that je"^ 7 ) (- iu) y l x l I = e-l cos( ^ 7)1 |aj|7 |x| is bounded for each 7 G R+ and thus 
J 7 (e~ s ^ (~ lw ) T l x l) is a tempered distribution and Lemma [4.21 can be applied. 
Using the same notion as in iflS"]). yields «*(— z) = s(7) (iz) 7 . 
Let z = re'* with r > and ip G (-tt,0). Then z G R + il" and Re((iz) 7 ) = 
cos(7(v? + tj)) |z| 7 , and thus the inequality in (B2) reads as follows 

(25) _ a ( 7 ) oos(7fa + |)) |z| 7 < C + TV log(l + |z|) (7 > 0) . 

This shows that condition (B2) is satisfied if and only if 5(7) cos(j((p + J)) > 0. 

a) For 7 G (0, 1] we get 5(7) = 1 and 7(1/3 + f ) G (-§, f ) for any ^ G (-7r,Q) 
and thus 5(7) cos(j(ip + ^)) > 0, i.e. condition (B2) is satisfied for any e > 0. 
Therefore .^(e^ 7 ) l x l) is causal for 7 G (0, 1]. 

b) Now we prove the non-causality of ^(e^ 7 ) ( " lw)T |x| ) for 7 G [3, 5]U[7, 9]U- • • . 
Since for these 7— values 5(7) = 1 we have to find a <p such that the sign 

of cos(7(<yS + |)) is negative. For 7 > 1 let < S < min T J, ^ and 

<Ps ■= (i-l) | + f Since W G (-f,-5) and 7 (w + f) = f +^ G (f, 3 ^), it 

follows that z := re' w G R + iK _ and cos(7(<^5 + ^)) < and hence condition 
(B2) cannot be satisfied for any e > 0. 

c) Now we prove the non-causality of ^(e^ 7 ) l x l) for 7 G (1, 3) U (5, 7) U 
(9, 11) U • • • . Since for these 7— values 5(7) = —1 we have to find a ip such that 
the sign of cos( 7 (<£>+-|)) is positive. For <p = — | it follows that z = re lip G it" 
(r > 0) and cos( 7 (<^+ ; |)) = 1, which implies at once that condition (B2) cannot 
be satisfied for any e > 0. This proves the Corollary. □ 

The following two Corollaries clarify for which values of 7 and luq the fre- 
quency power law a(u>) = «o M 7 together with the phase speed determined by 
the Kramers-Kronig relations ([!]), yield an attenuated spherical wave that sat- 
isfies the causality requirement. The derivation of the phase speed c = c(u>) 
corresponding to frequency power law with various 7— values can be found 

in na m mi m rm. 
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Corollary 4.4. Let ao and ujq be positive constants, 7 G M + \N and let the 
attenuation law a and the phase speed c be defined as follows: 

a(u>) = a |w| 7 , 

(26) 1 



c(w) c(uj ) 
Then a* defined as in (0) reads as follows 



-J- = a tan g 7 ) (M^-lwor 1 ). 



(27) a»(a;) = cos(f 7 ) +* a otan^-7j \u \ 



7-1, 



and £/ie wave defined by (QJ) /ias finite front speed t/76 (0, 1) and = 0. For 
t/iese cases i/ie wave /roni speed is equal to c(0). //^o = and 76(1, oo)\N or 
t^o > and 7 S R + \N, then the wave defined by (Q]) cannot have a finite front 



Proof. Let x e I 3 be arbitrary but fixed. From (|26|) and ([2|) we get at once 
identity ([27). 

i) Let wo = 0, then according to Lemma [4.31 Tie~ a * ' x ')(0 with a* defined as 
in (|27[) is causal if 7 € (0, 1) and non-causal if 7 S (1, oo)\N. 

ii) Now let uiq > 0. Then it follows that 



-Re(a*(— z)) = —01(7) |z| 7 +02(7) |Im(z)| for all z 6 

cos(f 7) 



where 01(7) := ao c ° 8( o ^t : ll , ^ G (— 7r, 0) is the argument of 2; and 02(7) 



a tan (§ 7) |w | 7 x . 

a) If 7 G (0, 1) then 01, 02 > and 

-Re (a*(— i 22)) = -01(7) Nl 7 + 02(7) (22] for all z 2 < , 

i.e. — Re (a*(— i 22)) growths like | ^2 i for sufficiently large — 2^. Since this term 
cannot be bounded by log(l + |z2 1)> condition (B2) in Lemma [4.21 cannot be 
satisfied for uj > and 7 G (0, 1). 

b) Now let 7 > 1 and 7 g" N. We note that 01(7) has the same sign as 
s(7) cos((tp + f)7), where 3(7) is defined as in Lemma l4~3l As in the proof 
of Lemma 14.31 one shows that for an appropriate choice of ip G (— tt, 0) the 
constant a± is negative, which shows that — Re (a* (— z(r))) growth like \z\ J for 
sufficiently large \z\. Therefore condition (B2) in Lemma l4~2l cannot be satisfied 
for ojq > and 7 G (0, 1). In summary we have shown that the front speed of 
the wave defined by (TTJ) and ([2"T|) cannot be finite if loq > and 7 G R + \N. 

iii) Now we show that the front speed for the case ojq = with 7 G (0, 1) is equal 
to c(0). We recall that if g a (x,t) is causal, then the front wave speed satisfies 
the condition v a < c{uiq). If the front wave speed v a at x is smaller than c(0) 
(ujq = 0), then 3 Q (x, t + S |x|) is causal for some S > 0. This means that there 
exist constants e > 0, C > and N > such that for all z G K + i (—00, — e): 

(28) -Re(a*(-z)) -Re(i(-z)5) < C + N log(l + |z|) . 

Since t^o = 0, we have 02 = 0. As above we get for all z — i Z2 G M~ 

-Re(a*(-z)) = - ai \z 2 p + 5\z 2 \ (ai > 0) , 
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which grows like \z2\ if <5 > and thus condition (f2"5| cannot be satisfied. 
This contradicts the fact that g a (x, t) is causal and hence we conclude that 
5 = 0. This proves that the wave front speed is equal to c(0) and concludes the 
proof. □ 



Now we come to the special case 7 = 1. 
Corollary 4.5. Let ao, luq > be constants and 

~a (-iw) 7 



a*(w) 



lim 

7—1 



(29) 
then 

(30) a(u>) = a \uj\ , 



cos(f 7) 



+ ia tan I — 7) |w | 7 1 u 



11 2 
— - - — — - = -a - log 



Ul 

w 



and the wave defined by ([!]) cannot have a finite front speed. Moroever, .?"(a*(co>)) 
is not causal. 



Proof. Definition 



implies 



a*(u>) = ao — iaoU) lim 



0+ cot(§ (1 + e)) ' 

Since both the numerator and the denominator of the last expression converge 
to zero, we can apply the rule of de I'Hospital, which yields 



a*(u>) — ao \ui\ = —iaoU) lim 



M e log M - ko| e log I w 1 



-o+ -f sin-(f (1 + e)) 
2 

= ia - us (log \w\ - log \oj \) . 

The result for the limit e — > 0— follows analogously. Comparing a* with ([2]) 
yields §U$ . For (x,y) e R 2 let 



u{x, y) := a Q \J x 2 + y 2 - a - y (log \J x 2 + y 2 - log \uj \) 

7T 



and 



v(x,y) := a - x (log \J x 2 + y 2 - log |w |) 



Then a*{z) = u(x,y) + iv(x,y) for z € C with a; = Re(z) and y = Im(z), and 
dU g^ y ^ dv< Qy yS> for every (x,y) G R 2 \(0,0). Since u(x,y) and v(x,y) do not 
satisfy the Cauchy-Riemann equations for every (x,y) € R 2 \(0,0), a*(— z) is 
not analytic on R + zR~, i.e. condition (Bl) in Lemma l4~2l is not satisfied. This 
shows that for fixed x e R 3 F(a*{u)){t) and J"( e - Q «^ l x l)(f) cannot be causal 
and concludes the proof. □ 
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5 Causality analysis of Szabo's wave equation 

The Green function of Szabo's equation (fTTj) for 7 > 0, 7 ^ N is given by (JTJ) 
with the following attenuation law and phase speed: 

1 \ is < < w a 1 1 Im(a*(w)) 

a(uj) = Rc(a*(w)) and — — = , 

c[lo) Co ijJ 

where 




(-t w) a + 2 a c - + t — 

cos (§7) c 



(cf. dispersion relation above equation pip .) Here the square root is understood 
as the primitive square root, since a(u>) has to be positive. Since for 7 G (0, 1) 



<**M » \ for|w|»l 

cos (f 7) 



and for 7 > 1, 7 £ N 



Q*H w ' for w|«l, 

cos (2 7) 

Szabo's model is a high frequency approximation of equation ([!]) if 7 G (0, 1) 
and as a small frequency approximation of equation ([9]) if 7 > 1, 7 ^ N. 

The next Proposition investigates the causality behaviour of Szabo's equa- 
tion. 

Proposition 5.1. Let 7 £ R + \N. TTie Green function G a of Szabo's equa- 
tion ill]) has finite front speed only 2/76 (0, 1). For 7 G (0, 1), the front speed 
of G a is cq. 

Proof. For the proof let x G R 3 be arbitrary but fixed. The Green function 
G(x, i) of wave equation (fTTj) has a front speed < Co < 00, if G(x,t + ^) is 
causal, i.e. if a*(w) defined by (|3"Tj) satisfies Lemma [4721 For convenience we set 
So := 



s (f t)' 



i) First we prove the Proposition for 7 G (0, 1). Since \J (i z) 2 + 2 6>o Co (i z) 7+1 
(primitive square root) maps R + iR _ analytically into R + iR~, a* (—2) maps 
R + iR~ analytically into R + iR". This proves condition (Bl). 
Let z = z\ + i Z2 with z\ G R, Z2 G R~ and 



! xi- 7 



B(z) := yi + 25 c J -1, 

then Re(a*(-2)) = Re(izB(z)) = -ziIm(B(z)) - z 2 Re(5(z)). Moreover, 
for Proposition 15.11 property (B2) in Lemma 14.21 reads as follows: there exist 
constants e > 0, C > and N > such that for all z G R + i (—00, — e): 

(32) — (*i Im (_B(z)) + z 2 Re (£(»)) < C + N log(l + |z|) . 

Co 
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We prove this inequality by showing that its left hand side is always negative 
or zero. 

a) First we prove that Re (B(z)) > if z 2 G R _ , which implies z 2 Re (B(z)) < 
if Z2 G M.~ . From 7 G (0, 1) and z 2 < 0, it follows that (i z) 1-7 has positive real 
part. Since the inversion of a complex number with positive real part yields a 
complex number with positive real part, 1 + 2 So Co 7 has a real part > 1. 
This implies that B(z) has positive real part. 

b) Now we show that z 1 lm(B(zj) < if z 2 G M~ and z\ G R. Let zi > 0, 
then the imaginary part of (iz) 1-7 is positive and since the inversion of a 
complex number with positive imaginary part yields a complex number with 
negative imaginary part, 1 + 2 So Co ( Jj) 7 has negative imaginary part. This 
implies that B{z) has negative imaginary part and therefore ziIm(_B(z)) < 
if Z2 G M~ and z\ > 0. Now let z\ < 0. Then the imaginary part of (i z) 1 1 is 
negative and since the inversion of a complex number with negative imaginary 
part yields a complex number with positive imaginary part, 1 + 2 So Co (j^) 7 
has positive imaginary part. This implies that B(z) has positive imaginary part 
and therefore zilm(B(z)) < if z 2 G M~~ and Z\ < 0. Clearly, if Z\ = then 
z\ Im (B(z)) — 0. In summary we have proven that the left hand side of (|32l) is 
smaller or equal to zero and thus the inequality holds. 

ii) The second part of the theorem is first proven for 7 € (1,3)\{2} and then 
for 7 > 3 with 7 ^ N. 

a) Let 7 G (1, 3)\{2}. (Indeed the following arguments hold for any 7 G (4n + 
l,4n + 3)\N with n S N .) Then a < and for z — i z 2 with z 2 < 0, 
condition ([32]) simplifies to 



Because 7 > 1 and So < 0, the term under the root is negative for sufncently 
large — z 2 and thus the real part of the root vanishes, which leads to the con- 
tradiction &± < C + N log(l + \z 2 \). Therefore condition (B2) cannot be valid 



b) Now let 7 > 3 and z — Z\ + i z 2 with z\ G M and z 2 < 0. We recall 
that B(z) = yjl + 2a c (iz) 7- - 1 - 1. Let z(r) = re iip with r > and 
(f := T^-j- — ^. Since 7 > 3, we have tp G (— f,0) and thus Zi > and 
z 2 < 0. Since (7 — 1) ((p + |) = 7r we have cos((7 — 1) (y? + §)) = —1 and 
sin ((7 - 1) (<p + f )) = 0. This shows that 

(33) 1 + 2 ao Co [i z(r)) 7 " 1 < for sufficiently large r. 

Therefore Re(_B(z)) = — 1 < 0, which together with z 2 < implies z 2 Re(_B(z)) = 
— z 2 > 0. This shows that the first left hand side term of ((3"2"f is positive. More- 
over, (|33jl implies that lm(B(z)) = lm(B(z) + 1) > for sufficiently large — z 2 . 
From this together with z\ > we obtain ii lm(B(z)) > 0, which shows that 
the second term on the left hand side of (f32j) is positive, too. If r is sufficiently 

-7 — 1 

large then Im(B(z(r))) is of the order of with 7 > 3 which cannot be 

bounded by C + N log(l + r). Hence condition (B2) cannot be true for 7 > 3. 
In summary we have shown that the Green function of Szabo's equation cannot 
be causal for 7 e (1, oo)\N. 

iii) Now we show that the front speed of G a is cq, if 7 e (0, 1). Since G(x, t+ 



— - - Re [Jzl + 2 5 c (-z 2 ) 7+1 < C + N log(l + |z 2 |) . 




for any 7 eTl,3)\{2}. 
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is causal, the front wave speed v a satisfies v a < cq. If v a at x is smaller than 
Co, then G(x, t + ^ + S |x|) is causal for some 5 > 0. This means that there 
exist constants e > 0, C > and N > such that for all z G K + i (-co, — e) 
(cf. flU in Corollary [ 



-Re(iz (P(z) - <5)) < C + N log(l + \z\) 



since 



He(o,(-«)) + Re(t (-«) 5) = -Re(iz - <5)) . 
For the setting z — i z 2 with Z2 < 0, we obtain 



Re |z 2 | 



<5 + 1 - 4/ 1 + 2 S c 



1-7 



< C + N log(l + |z 2 |). 



which cannot be true if \z2\ is sufficiently large. Hence we conclude that 5 = 0. 
This proves that the wave front speed is Cq and concludes the proof. □ 



6 Causality analysis of the thermo- viscous wave 
equation 

The operator of the thermo-viscous wave equation (cf. e.g. [7]) 

has order 3 and principal part P3 (X) = To t ^2j=i x ji where X := (t, x±, x<z, xs) T £ 
R 4 . Since P 3 (N) = for N := (1,0,0,0) T , the plane {X e M 4 | (X, N) = 0} 
in M 4 is characteristic with respect to P(D). According to Theorem 8.6.7 in |5j 
the thermo-viscous wave equation with vanishing source term has a solution 
Ptv € C°°(IR 4 ) such that supp(p tlI ) = M 3 x Kq . This shows that the Green 
function of the thermo-viscous wave equation is not uniquely determined. (The 
existence is guaranteed, since P(X) := (1 + r t) y^ =1 — ^ is not the zero 
polynomial.) Theorem 16.21 below shows that the Green function of the thermo- 
viscous wave equation cannot have a finite front speed and that a solution of 
the thermo-viscous wave equation depends on its history. This explains why its 
Cauchy problem has no unique solution. For this theorem and Theorem 1 7. II we 
need the following lemma. 

Lemma 6.1. Let S + = {/ e S(R) | supp(f) C R+} and 7 G (1,2]. The time- 

1 

convolution operator T 7 2 defined by the kernel 

k i(t) :=^LW 1 = ) (t) 

is an isomorphism ofS+. Here the square root is understood as the primitive 

1 x 
square root. The inverse of T 7 2 is the time-convolution operator Lj with the 

kernel 

k 1 

Again the square root is understood as the primitive square root. 



■ (<) := 4= T + (-irouF-A (t) . 

V 2 7T v ' 
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Proof. First we show that the kernel of T 7 and L 2 have supports in [0, oo). 
Let B(z) := yfl + (jto-z) 7-1 andzeM + sR". Then 



supp Ife ij C [0,oo) and supp Ifc ij C [0, 00), 

if and only if u\{z) :— B(z) and U2{z) :— -g^j satisfy Theorem 14.11 Since 

v/1 + (itq z)" 1-1 maps R + i R~ analytically into [1, 00) -MR and 1 + (ito z) 7_1 
cannot vanish on R + iK _ , Mi(z) and 112(2) map R + iR~ analytically into 
[1, 00) + i M. Therefore property (Al) in Lemma l4~2l is satisfied. We have for all 
z G R + iR~ with \z\ » 1: 



|S(z)| < Ci |z| (7 ~ 1)/2 and 



B(z) 



1 \ (7-l)/2 



<Ci (7 6(1,2]) 



for some constants C\, C2 > 0. Therefore property (A2) in Theorem 14.11 is 

j_ 1 

satisfied for ui(z) and 1*2(2), which proves that the kernels of T 7 and L 7 have 
support in [0, 00). 
Since 

T 7 2 f G T>' + for every / G S+ , 

and 



1 (t| /) ( W ) = /M = G 5 for every / G <S+ , 

V / ^1 + (-irow)T- 1 



the convolution operator T 7 is well-defined on 5+. Since 5 is invariant under 

1 

multiplication by a polynomial, it follows analogously that L 2 maps S+ into <S+ 
and is well-defined. Since the Fourier transform is an isomorphism on S and 

t- 1 (4 4 /) =f = T~ 1 (t| 4 /) fOT cver y / e 5 . 

it follows that L 2 : 5+ — * 5+ is the inverse of T 7 2 : X + — > 5+ and T 7 is an 
isomorphism of <S+. □ 

With the help of the Laplace transform table in [5j (cf. Appendix 2), the 

1 

kernel of T 2 2 can be calculated as follows: 

kr (t) = -jL T ( . \ ) (t) = V^C- 1 ( n 1 ) (t) 



— e-*H{t). 

T t 

where H (t) denotes the Heaviside function. This shows that for 7 = 2 the kernel 

1 1 1 

of T 2 2 decreases exponentially. In the following we denote T 7 T 7 by T 7 . 



Proposition 6.2. Let to, Co > be constants and A(u>) := 1 + y/l+ljo oj) 2 
for any 10 G R. The Green function of the thermo-viscous wave equation 

9\ n2 ,,,, l9 2 G(x,i) sf . sf . m3 



(34) 1 J + T - 1 V 2 G(x, t) - ? ^2 ' y = -8(x) 5(t) x 



1G 



cannot have a finite front speed and is given by G = T 2 p a , where p a is defined 
as in {TJj with attenuation law 

(35) a(w) 



y/2A(u) (A(ui) - 1) c 
and phase speed 

(36) c(w) = == c . 

y/A(u) 

Proof. Applying the inverse Fourier transform to the thermo- viscous wave equa- 
tion yields 

V 2 G(x,w) + k 2 (uj)G(x,uj) = S ^ with 

This problem has the solution G(x, uj) = , 1 r e . , , . where the square 

1 V^jrtl-iTju) 47r|x| ' 

root of 1 — i To a; is understood as the root with positive real part. We assume 
that G satisfies the causality requirement (j!4|) which in particular implies that 
supp(G(x, ■)) £ [0, oo). Then the Green function can be written as follows 

1 f e -i{<*it-K u ) l x l) 

( :is ) G(x, t)=T 2 | — J _-j duj | =: T 2 p Q (x, t) 

where T2 denotes the time-convolution operator in Lemma 16. II for 7 = 2. The 
last identity is equivalent to 

Id + t — J G(x, i) = p Q (x, t) , 

which shows that p a has finite front speed if and only if G has finite front speed. 
Let c\ > be arbitrary but fixed. We prove a contradiction by showing that 
p a (x, t + cannot be a causal distribution for any x £ M 3 . Comparing l|38p 

and ([T]) shows that a(u>) — Jm(k(tj)) and = -^ e W") , Since a > is 
required, the imaginary part of k(cu) must be positive and therefore we choose 
the positive sign for fc, i.e. 



, / % W 1 / A(u>) TqUJ 

(39) fc( li ,) = A /_Li + i . 



c — 1 V V 2 v / 2A(lu) 



with A(w) := 1 + ^/l + (t cj) 2 . From this we get the attenuation law (|35]) and 
the phase speed (131)1) . If p Q (x, t+^) is causal, then property (B2) in Lemma HT!? 
must be satisfied for —i k(uj) + i i.e. for each xel 3 there exist constants 
e > 0, G > and N > such that for all z £ R + 1 (-00, -e): 

Re(ifc(-z))+Re N-^-J < G + TV log(l + \z\) . 
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For z — i z 2 with 22 < we get 



1 12 <C + N\og(l + \z 2 \). 



VI - T Z 2 C Cl 

which cannot be true for sufficiently large — z 2 and finite c\ > 0. Hence property 
(B2) cannot be satisfied for a* of the wave p a if c\ is finite. This contradiction 
proves that the front speed of G cannot be finite and concludes the proof. □ 



7 Causal wave equations obeying attenuation power 
laws with 7 G (1,2] for small frequencies 

As we have seen in Section [4j the frequency power law for 7 > 1 cannot hold. 
Now we present admissible attenuation laws that permit approximate frequency 
power laws with 7 G (1,2] for sufficiently small frequencies. 

For given constants 76 (1,2], 0<ci <oo and r > we define 

(40) M == 7t==^| w _ f ' 

ci V 1 + H T o w ) 7 

where the square root is understood as the primitive square root. This implies 
for the attenuation law: 

( \ ~ \ 17 ^ Sin (f (7-1)) 

ctyuj) ~ ao bo w with «o = 

T Cl 
1 1 

for sufficiently small frequencies. Moreover, let T 7 2 and L 7 be defined as in 
Lemma I(TT1 and let the operators -D* and -D*,|x| be defined as in Section [3] Then 

1 id 

D* = —T*— and D, jW = 

Ci OT 11 

and wave equation (|22[) (with vb replaced by cq) reads as follows 



(41) W 



1 

Cl 



i 2 1 ^2 



Id + — T 7 2 



1 5V 
9i 2 



For 7 = 1 we obtain the classical wave equation without damping and for 
7 = 2, we obtain a modified thermo- viscous wave equation. Since Re(a*) is 
equal to (|35|) for 7 = 2, the modified thermo- viscous wave equation obeys for 
7 = 2 the same attenuation law as the thermo-viscous wave equation (if the 
source term / is replaced by L 2 /). 

Proposition 7.1. The Green function of wave equation has finite and 
constant wave front speed cq . 

Proof. For the proof let x € M 3 be arbitrary but fixed. The Green function 
G(x, t) of wave equation (|4T|) has a front speed < Co < 00, if Gt v (x,t + ^) is 
causal, i.e. if a*(u>) defined by (HP)) satisfies Lemma [4.21 We recall that the 
square root in the definition of a* is understood as the primitive square root. 
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According to the proof of Lemma f6.1l a * ~ z satisfies property (Bl) and thus 
a*(— z) satisfies property (Bl), too. 

Property (B2) in Lemma [4.21 reads as follows: there exist constants e > 0, 
C > and N > such that for all z G R + i (— oo, -e): 

(42) -Re(a.(-«)) = -Re( ' Z ) < C + N log(l + |z|) . 

Therefore (B 2) is satisfied if -Rje(a*(-z)) < for each z G R + iR _ . Let 
B(z) := ci + (iToz)''- 1 and z = zi + i z 2 G R + t R~ then 

-Re (a* (-.?)) = z 1 Im(B- 1 (z)) + z 2 Re {B~ 1 {z)) . 

Since 7 G (1,2] and z 2 < 0, it follows that (jroz) 7-1 has positive real part 
and thus B(z) has also positive real part. The inversion of a complex number 
with positive real part yields a complex number with positive real part and thus 
Re {B- 1 {z)) > for any z 2 < 0. This proves that 

(43) z 2 Re {B- 1 {z)) < for any z G R + i Mr . 

For z\ > the imaginary part of {itqz) 1 ' 1 is positive and thus lm(B(z)) > 
0. The inversion of a complex number with positive imaginary part yields 
a complex number with negative imaginary part and thus Im(.B -1 (z)) < 0. 
Therefore we infer that z x Im (B~ 1 (z)) < for any z G R + + i R~ with z\ > 0. 
For z\ < the imaginary part of (zto z) 7 " 1 is negative and thus lm(B(z)) < 0. 
Since the inversion of a complex number with negative imaginary part yields a 
complex number with positive imaginary part we conclude that Im (_B _1 (z)) > 
0. Hence z\ Im ( y B^ 1 (z)) < for any z G R~ + i R~ with z x < 0. For z\ = we 
get z\ Im (_B _1 (z)) = 0. In summary we have proven that 

(44) ziIm(S" 1 (z)) < for any z G R + i R~ . 

(|4"3"|) and (|4^|) imply that the left hand side of (|4"2"]) is always non-positive and 
therefore (05]) is true. This shows that G(x, t) has a front speed < c < 00. 

Now we show that the front speed of G is equal to cq. If the front speed is 
v a (x) < Co for any x G R 3 , then (|42|) must hold for 

a*(— z) := a*(— z) + ie (— z) . 

For z := i z 2 with sufficiently large — z 2 we obtain 

-Re (a* (-z)) = = = + e (-z 2 ) , 

ci Vl + (-^o^) 7 - 1 

which is positive and of the order This shows that condition (B2) can only 
be true if e — 0. This concludes the proof. □ 

Remark 7.2. For j = 2 let G Co denote the solution of wave equation J^i[ ) anc? 
let Gtv denote the solution of the thermo-viscous wave equation {3$ with Co 
replaced by c\. Then one can show that 

lim L G CQ (x, t) = G tv (x, i) for each x G R 3 and t G R + , 

Co — >oo 

which shows again that the front speed of Gtv is infinite. 



19 



8 Acknowledgement 



This paper was partly supported by the " Frankrcichschwerpunkt" of the Univer- 
sity of Innsbruck. I would like to thank Xavier Bonnefond and Pierre Marcchal 
for fruitful discussions. 



References 

E.J. Beltrami and M. R. Wohlers. Distributions and the Boundary Values 
of Analytic Functions. Academic Press, New York and London, 1966. 

Leon Brillouin. Wave propagation and group velocity, volume 3 edition. 
Academic Press, New York and London, 1960. 

R. P. Feynman, R. B. Leighton, and M. Sands. The Feynman Lectures on 
Physics. Vol 1-2. Addison- Wesley Publishing Company, Reading, 1963. 

Harro Heuser. Gewohnliche Differentialgleichungen, volume 3 Auflage. 
B.C. Tcubncr, Stuttgart, 1995. 

L Hormandcr. The Analysis of Linear Partial Differential Operators I, 
volume 2nd Edition. Springer- Verlag, New York, 2003. 

A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo. Theory and Applications 
of Fractional Differential Equation, volume Mathematics Studies 2004. El- 
sevier, New York, 2006. 

L. E. Kinsler, A. R. Frey, A. B. Coppens, and J. V. Sanders. Fundamentals 
of Acoustics. John Wiley and Sons, Inc., New York, 2000. 

S. Lang. Complex analysis. Springer- Verlag, New York, 1993. 

I. Podlubny. Fractional Differential Equations, volume Mathematics in 
Science ans Engineering 198. Academic Press, New York, 1999. 

T. D. (Ed.) Rossing. Springer Handbook of Acoustics. Springer Verlag, 
New York, 2007. 

T. L. Szabo. Time domain wave equations for lossy media obeying a fre- 
quncy power law. J. Aucoust. Soc. Am., 96(1), 1994. 

T. L. Szabo. Causal theories and data for acoustic attenuation obeying a 
frequncy power law. J. Acoust. Soc. Am., 97(1), 1995. 

K. R. Waters. On a time-domain representation of the Kramers-Kronig 
dispersion relation. J. Aucoust. Soc. Am., 108(5), 2000. 

K. R. Waters, M. S. Hughes, G. H. Brandcnburgcr, and J. G. Miller. On the 
applicability of Kramers-Kronig relations for ultrasonic attenuation obey- 
ing a frequaency power law. J. Aucoust. Soc. Am., 108(2), 2000. 

K. R. Waters, M. S. Hughes, J. Mobley, and J. G. Miller. Differential 
forms of the Kramers-Kronig dispersion relations. IEE Trans. Ultrason., 
Ferroelect., Freq. Contr., 50(1), 2003. 



20 



[16] K. R. Waters, J. Mobley, and J. G. Miller. Causality-imposed (Kramers- 
Kronig) relationships between attenuation and dispersion. IEE Trans. Ul- 
trason., Ferroelect.,Freq. Contr., 52(5), 2005. 



21 



